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ABSTRACT 

Context. The Bulge is the least understood major stellar population of the Milky Way. Most of what we know about the formation 
and evolution of the Bulge comes from bright giant stars. The underlying assumption that giants represent all the stars, and accurately 
trace the chemical evolution of a stellar population, is under debate. In particular, recent observations of a few microlensed dwarf stars 
give a very different picture of the evolution of the Bulge from that given by the giant stars. 

Aims. We aim to resolve the apparent discrepancy between Bulge metallicity distributions derived from microlensed dwarf stars and 
giant stars. Additionally, we aim to put observational constraints on the elemental abundance trends and chemical evolution of the 
Bulge. 

Methods. We perform a detailed elemental abundance analysis of dwarf stars in the Galactic bulge, based on high-resolution spectra 
that were obtained while the stars were optically magnified during gravitational microlensing events. The analysis method is the same 
as for a large sample of F and G dwarf stars in the Solar neighbourhood, enabling a fully differential comparison between the Bulge 
and the local stellar populations in the Galactic disc. 

Results. We present detailed elemental abundances and stellar ages for six new dwarf stars in the Galactic bulge. Combining these 
with previous events, here re-analysed with the same methods, we study a homogeneous sample of 15 stars, which constitute the 
largest sample to date of microlensed dwarf stars in the Galactic bulge. We find that the stars span the full range of metallicities from 
[Fe/H] = -0.72 to +0.54, and an average metallicity of ([Fe/H]} = -0.08 ± 0.47, close to the average metallicity based on giant 
stars in the Bulge. Furthermore, the stars follow well-defined abundance trends, that for [Fe/H] < are very similar to those of the 
local Galactic thick disc. This suggests that the Bulge and the thick disc have had, at least partially, comparable chemical histories. 
At sub-solar metallicities we find the Bulge dwarf stars to have consistently old ages, while at super-solar metallicities we find a wide 
range of ages. Using the new age and abundance results from the microlensed dwarf stars we investigate possible formation scenarios 
for the Bulge. 

Key words, gravitational lensing — Galaxy: bulge — Galaxy: formation — Galaxy: evolution — stars: abundances 



1. Introduction 

The Galactic bulge is a major stellar component of the Milky 
Way. Esti mations of its m ass range from 10% (lOorj 119771) 
to 25 % dSofue et all 120091) of the total stellar mass of the 
Galaxy. It is a peanut shaped barred bulge and occupies the in- 
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ner ~1 kpc of the Galaxy (lFrogellll988l) . Recent results for the 
Bulg e shape give a scale length for the bar major axis of ~ 1 .5 kpc 
(Rattenbury et al. 2007). The kinematic properties of the Bulge 
are intermediate between a rotationally supp orted system and a 
veloc ity dispersion dominated system (e.g., iMinniti & ZoccalJ 
2008). The markedly different stellar populations that inhabit 
this region of the Galaxy make it important to discern the for- 
mation and evolution as part of the understanding of the overall 
formation of the Galaxy. In addition, bulges are common fea- 
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tures among luminous galaxies, and determining the origin of 
the Bulge and the corresponding observational signatures (e.g., 
boxy isophotes, high [a/Fe] ratios, no age dispersion) are impor- 
tant steps in decoding the formation of bulges in general. 

The differences among the Bulge, the thin and thick discs 
in the Solar neighbourhood, and the Galactic stellar halo have 
been discovered and detailed through extensive photometric and 
spectroscopic observations of stars. Red giants are generally the 
only stars bright enough for high-resolution spectroscopy at the 
distance of the Bulge, and they have been intensively studied 
in the optical and infr ared, particularly in a few fields including 
Baade's window (e.g. McWilliam & Rich 1994; Fulbright et al. 
2006: lZoccali et al ll2008tlMelendez et al.l2008tlcunha & S mith 
2006: iRvde et al l 120091) . Chemically, the Bulge is decidedly 



more metal-rich t h an the stellar halo, with the mean [Fe/H] from 
the IZoccali et al.l (|2008) sample slightly below solar metallic - 
ity, and a possible vertical metallicity gradient of 0.6 dex per 
kpc, although a comparison with M giants in the inner Bulge 
may indicate a fl attening of the gradient in the central regions 
dRich et al.H2007h . The abundance ratios show enhanced [a/Fe] 
ratios that persist to higher [Fe/H] than in loca l thin disc stars, 
but show good agreement with thick disc giants ( Mel endez et al.1 
2008; lAlves-Brito et al.l 120101) . Furthermore, the stars at the 
Bulge main-sequence turnoff are red, and isochrone fitting has 
shown that the major i ty of the stars in the Bulge are old (e.g., 
Holtzman et all 119931; iFeltzing & Gilmord 120001 IZoccali et ail 
20031: IClarkson et al.ll2008h ~ 



These observations have fuelled an intense debate on the ori- 
gin of the Bulge. The established concordance model of cosmol- 
ogy provides the framework to understand the formation of the 
Milky Way by hierarchical merging, but the formation of the 
components of the Galaxy requires important physics, such as 
star formation and feedback, that is unresolved and parametrised 
in current models. However, there are two basic scenarios by 
which a bulge is built up in simulations (see lRahimi et al.l [2009. 
and references therein). The first is by mergers, where sub- 
clumps coming together in the early phases of Galactic evolu- 
tion combine to form the Bulge out of both accreted stars and 
stars formed in situ. The second is secular evolution, where 
the Bulge is created gradual ly out of the Galactic disc (e.g., 
iKormendv & Kennicutti2 004). The merger model is favoured by 
the metallicity gradient, while the secular evolution model is fa- 
vored by cylindrical rotation and by agreement (in terms of mean 
metallicity abu ndance trends, and ages) with the Galactic thick 
disc (e.g jHowardet al.ll2009l) . 

As mentioned above, spectroscopic observations of stars in 
the Bulge have usually been confined to giants. This limits our 
knowledge of the Bulge in several ways. First, because much 
of our kn owledge of the Solar ne ighbourhood relies on dwarf 
stars (e.g. Edvardsson et al. 1993), any systematic offsets be- 
tween the metallicity scal e of giant and dw a rf star s is cause for 
concern. In this context, iTavlor & Croxalll (120051) have shown 
that there is a lack of very metal -rich ([Fe/H] > +0.2) gi- 
ants in the Solar neighbourhood and that the mean metallic- 
ity of local gi ants is lower than for dwarfs. Recent s tudies of 
nearby giants dLuck & Heiter 2007; Takeda et al. 2008) confirm 
the lack of very metal-rich stars. ISantos et al.ld2009l) suggest that 
there may be systematic errors in the metallicity determinations 
of metal-rich giants. Secondly, there are several pieces of evi- 
dence that giant stars may not accurate ly represent all the stars. 
For the metal-rich cluster NGC 679 1 . iKalirai et all (120071) pro- 
posed that the explanation for the large number of low-mass 
He white dwarfs was that about 40 % of the stars do not be- 
come red clump stars, skipping entire phases of stellar evolu- 



tion altogether because of high mass loss. lKilic et al.l (f2007) pro- 
posed a similar mechanism to explain the non-binary He white 
dwarfs found in the Solar neighbourhood. Finally, with turnoff 
stars, ages can be determined for individual stars, and an age- 
metallicity relationship derived, which is not possible for gi- 
ant stars. This means that a true study of the Bulge requires 
the study of dwarf stars. However, at the distance of the Bulge, 
dwarf stars are too faint for abundance analyses based on high- 
resolution spectra. Turnoff stars in the Bulge have V magnitudes 
ar ound 19 to 20 (compare, e .g., the colour-magnitude diagrams 
in Feltzing & Gilmore 2000). However, in the event that a Bulge 
dwarf star is lensed by a foreground object, the magnitude of the 
star can increase by more than 5 magnitudes, in which case a 
high-resolution spectrum can be obtained and the star analysed 
in a similar manner as the dwarf stars in the Solar neighbour- 
hood. 



There were several spectroscopic observations of mi- 
crolensed Bulge stars in the 1990s, but th e first high-resolutio n 
spectrum of a dwarf star was presented in Min niti et al.l (Q998). 
Complete high-resolution spectroscopic abundance analyses 
have b een published for eight mic r olensing events toward the 



Bulge (|Johnson et al 



2007, 2008 



Cohen et all 120081 [2009b; 



iBensbv et al.ll2009albi: Epstein et all2009l) . Initially, it appeared 
that the microlensed Bulge dwarf stars were ty pically much 
more metal-rich than the giants in the Bulge, and Epstein et al. 
(2009) found, using a Kolmogorov-Smirnow (KS) test, a very 
low probability of only 1.6% that these eight microlensed dwarf 
stars in the Bulge were drawn from the same met allicity distri- 
bution (MDF) as the sampl e of Bulge giants from IZoccali et al.l 
(20081). ICohen et aTl (120081) proposed that a similar mechanism 
to the one occurring in NGC 6791 was occurring in the Bulge, 
preventing the most metal-rich stars from being included in 
the giant sur v eys. Arguments against this idea presented by 
IZoccali et all {2008) were based on the luminosity function 
along the main sequence and red giant branch, showing no lack 
of RGB stars, with respect to the prediction of theoretical mod- 
els. 



In addition, the microlensed Bulge dwarf stars showed good 
agreement in abundanc e ratios with the thi ck disc stars in the 
Solar neighbourhood (Bens bvet al.l l2009allbl) . However, com- 
parisons have been hampered because of the small number of 
microlensed stars that did not always cover the [Fe/H] range 
of interest. Also, individual age estimates were provided for 
stars near the turnoff subgiant branch, including some stars 
that could be younger than the canonical old Bulge population 
dJohnson et al.ll20oilBensbv et al.ll2009bl) . 



We will 
suits for 
bulge. We al so 

Caval lo et al l d2003l) : ICohen et al.l (|2009b) and 
(2009). Comb ining these data with the results from 
(2009a) and Bensb y et al. (l2009bl) (which includes a re-analysis 



here present detailed elemental abundance re- 
six new microlensing events toward the Galactic 
re-an alyse the ev e nts pre viousl y studied by 



Epstein et al 



Bensbv et al 



of the events from lJohnson et al.l2007l 120081: ICohen et al.ll2008l) 
we now have a sample of 15 microlensed dwarf stars in the Bulge 
that have been homogeneously analysed using the exact same 
methods. 
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Table 1. Summary 1 of the, so far, 15 dwarf star microlensing events in the Bulge that have been observed with high-resolution 
spectrographs. They have been sorted according to their metallicities (as given in Table[2]). 



Object 


RAJ2000 


DEJ2000 


I 


b 


T E 


T max 


Afnax 


Tobs 


Exp. 


S/N 


Spec. 


R 


Ref. 




[hh:mm:ss] 


[dd:mm:ss] 


[deg] 


[deg] 


[days] 


[HJD] 




[MJD] 


[s] 










OGLE-2009-BLG-076S 


17:58:31.9 


-29:12:17.8 


+ 1.21 


-2.56 


36.9 


4916.46 


70 


4916.291 


7200 


30 


UVES 


45 000 


B09a 


MOA-2009-BLG-493S 


17:55:46.0 


-28:48:25.8 


+ 1.25 


-1.84 


13.2 


5094.61 


150 


5093.980 


7200 


40 


UVES 


45 000 


TW 


MOA-2009-BLG-133S 


18:06:32.8 


-31:30:10.7 


+0.05 


-5.19 


26.6 


4932.10 


70 


4932.233 


7200 


35 


UVES 


45 000 


TW 


MOA-2009-BLG-475S 


18:02:27.3 


-27:26:49.9 


+3.16 


-2.45 


34.0 


5084.92 


62 


5084.980 


7200 


20 


UVES 


45 000 


TW 


MACHO- 1999-BLG-022S* 


18:05:05.8 


-28:34:39.5 


+2.46 


-3.51 


265.0 


1365.50 


28 


1366.315 


12600 


80 


HIRES 


29 000 


C03 


OGLE-2008-BLG-209S 


18:04:50.0 


-29:42:35.3 


+ 1.44 


-4.01 


19.5 


4606.09 


30 


4606.833 


5400 


30 


MIKE 


47 000 


B09b 


MOA-2009-BLG-489S 


17:57:46.5 


-28:38:57.8 


+ 1.61 


-2.14 


58.9 


5095.52 


103 


5094.982 


7200 


65 


UVES 


45 000 


TW 


MOA-2009-BLG-456S 


17:48:56.3 


-34:13:32.3 


-4.16 


-3.34 


36.2 


5090.94 


77 


5090.982 


7200 


45 


UVES 


45 000 


TW 


OGLE-2007-BLG-514S 


17:58:03.0 


-27:31:08.2 


+2.62 


-1.63 


76.0 


4386.28 


1200 


4385.985 


4800 


30 


MIKE 


25 000 


E09 


MOA-2009-BLG-259S 


17:57:57.6 


-29:11:39.1 


+ 1.15 


-2.45 


69.1 


5016.77 


220 


5016.227 


7920 


50 


UVES 


45 000 


TW 


MOA-2008-BLG-311S 


17:56:53.7 


-31:23:40.3 


-0.87 


-3.35 


15.7 


4655.40 


400 


4654.952 


7200 


85 


MIKE 


29 000 


C09 


MOA-2008-BLG-310S 


17:54:14.4 


-34:46:37.7 


-4.09 


-4.56 


7.1 


4656.39 


280 


4655.957 


7200 


90 


MIKE 


41000 


C09 


OGLE-2007-BLG-349S 


18:05:23.0 


-26:25:27.1 


+4.38 


-2.52 


109.4 


4348.56 


400 


4348.237 


3050 


90 


HIRES 


48 000 


C08 


MOA-2006-BLG-099S 


17:54:10.2 


-35:13:34.9 


-4.48 


-4.78 


30.1 


3940.35 


380 


3940.090 


2400 


30 


MIKE 


19 000 


JOS 


OGLE-2006-BLG-265S 


18:07:18.9 


-27:47:44.0 


+3.38 


-3.55 


28.6 


3893.24 


212 


3892.581 


900 


45 


HIRES 


45 000 


J07 



* Given for each microlensing event is: RA and DE coordinates (J2000) read from the fits headers of the spectra (the direction where the telescope pointed during observation); galactic 
coordinates (/ and b)\ duration of the event in days (Te); time when maximum magnification occured (T mx ); maximum magnification (AmaxY* ume when event was observed with 
high-resolution spectrograph (T i, s ); magnification at the time of observation (A />, f ); the exposure time (Exp.) , the measu red signal-to-noise r ation per pixel at ~6 400 A; the spectrogr aph 
that w as used; the spe ctral resolut ion; a nd the reference whe re the star first appeared: TW=This Work. B09a= Bensby et al. 1 2009a). B09b= Bensby et al. i 2009b), J07= Johnso n et alj 
<2007l) . JO MJohnson etai] UtM) . C08= ICohen etal] (2003) , C09= ICohen etall I2009bl) . C03= )Cavallo etai](2003l) . 

* MACHO ID: 109.20893.3423. 



2. Observations and data reduction 

In order to trigger observations of these highly magnified stars 
we rely on the OGLEQ and MOA0 projects that every night mon- 
itor about 100 million stars toward the Bulge to detect variations 
in their brightnesses. If an object shows a well-defined rise in 
brightness, a microlensing alert is announced. Every year, 600- 
800 events are detected. Based on the photometric data obtained 
by the MOA and OGLE surveys it is possible to model the event, 
and make predictions of the length of the event, peak brightness, 
and time of peak brightness. Stars are identified as likely dwarf 
stars based on their unlensed magnitudes and colour differences 
relative to the red clump stars. This is done in instrumental mag- 
nitudes. The majority are low magnification events, and only a 
few have unlensed brightnesses of V — 18 - 20, characteristic of 
dwarf stars in the Bulge at a distance of ~ 8 kpc. During a reg- 
ular Bulge season in the Southern hemisphere, typically around 
10 high-magnification events of dwarf stars in the Bulge are de- 
tectecQ. To catch these unpredictable events, we have an ongo- 
ing Target of Opportunity (ToO) program at the ESO Very Large 
Telescope on Paranal in Chile. Observations can then be trig- 
gered with only a few hours notice. 

2.1. The 2009 events observed with UVES 

On March 21, 2009, the OGLE early warning system identified 
OGLE-2009-BLG-076 to be a possible high-magnification mi- 
crolensing event toward the Galactic bulge. As the magnitude of 
the source star before the microlensing event indicated that it was 

1 OGLE is short for Optical Gravitationa l Lens Experiment, 
|http : //ogle . astrouw . edu . pi |dUdalski|[2003h . 

2 MOA is short for Microlensing Observations for A strophysics, 
http://www.phys.canterbury.ac.nz/moa (e.g., Bon d et ail 
1200 lh . 

3 With the new OGLE-IV camera that will be in operation in the first 
half of 2010 the field-of-view will increase from the 0.3 square degrees 
of OGLE-III to 1.4 square degrees, resulting in a substantial increase in 
the number of detected microlensing events. 



either a dwarf or a subgi ant star, we tri ggered observations with 
the UVES spectrograph (iDekker et al.ll2.Q00l) located on UT2 at 
the ESO Very Large Telescope on Paranal. Due to the limited 
visibility of the Bulge in March, the target had to be observed to- 
wards the end of the night. Hence, OGLE-2009-BLG-076S was 
observed on March 26, a few hours after reaching peak bright- 
ness (see Fig. [T}. A few weeks later, on April 11, we observed 
the MOA-2009-BLG- 133, this time first alerted by the MOA col- 
laboration, also with UVES. The third event was observed in the 
beginning of July, MOA-2009-BLG-259S. This object was ob- 
served during the UVES red arm upgrade, so we could only ob- 
tain a spectrum with the blue CCD that has a limited wavelength 
coverage of 3700-5000 A. In September, the end of ESO ob- 
serving period P83, we saw an explosion of microlensing events 
and another four source stars were observed with UVES: MOA- 
2009-BLG-475S on Sep 10, MOA-2009-BLG-456S on Sep 16, 
MOA-2009-BLG-493S on Sep 19, and MOA-2009-BLG-489S 
on Sep 20. The UVES red arm was now back with two new red 
CCDs with increased sensitivities. 

For the 2009 observations with UVES listed above, each tar- 
get was observed for a total of two hours, split into either four 30 
min or three 40 min exposures. Using UVES with dichroic num- 
ber 2, each observation resulted in spectra with wavelength cov- 
erage between 3760-4980 A (blue CCD), 5680-7500 A (lower 
red CCD), and 7660-9460 A (upper red CCD). In all cases a slit 
width of 1" was used, giving a resolving power of R « 45 000. 

On all occasions, right before or right after the observations 
of the microlensed targets in the Bulge, a rapidly rotating B star, 
either HR6141 or HR8431, was observed at an airmass similar 
to that of the Bulge stars. The featureless spectra from these B 
stars were used to divide out telluric lines in the spectra of the 
Bulge stars. Also, at the beginning of the night of April 1 1 we 
obtained a solar spectrum by observing the asteroid Pallas. 

Data taken before the upgrade of the UVES red CCD in mid- 
July were reduced with the UVES pipeline (CPL version 3.9.0), 
while the data taken after the upgrade were reduced with version 
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Fig. 1. Light curves for the eight new microlensing events. The photometry comes from the surveys indicated by their names (MOA 
or OGLE), except for MACHO- 1999-BLG-022S that has data from from both OGLE (circles) and binned MACHO data (crosses). 
Each plot has a zoom window, showing the time intervals when the source stars were observed with high-resolution spectrographs. 
In each plot the un-lensed magnitude of the source star is also given (Is). 



4.4.5. Typical signal-to-noise ratios per pixel at 6400 A are given 
inTableQ] 

The light curves for the seve n microlensing even ts (includ- 
ing OGLE-2009-BLG-076S from Ben sby et al.f2009al) observed 
with UVES in 2009 are shown in Fig. [1] in which we have also 
indicated the time interval during which they were observed with 
high-resolution spectrographs. Positions, amplifications, times 
of observation, and exposure times are given in Table Q] 

2.2. Los MACHOs 

Cavallo et alj d2003h presented the first detailed elemental abun- 
dance study of microlensed dwarf stars in the Bulge. Their 
analysis was of a "preliminary" nature, so we decided to re- 
analyse the stars that they classified as either dwarf or sub- 
giant stars. There are four such stars: MACHO- 1997-BLG- 
045S, MACHO-1998-BLG-006S, MACHO-1999-BLG-001S, 
and MACHO- 1999-BLG-022S. 

The observations of these stars were carried out from 1997 
to 1999 with the HIRES spectrograph on the Keck I telescope on 
Hawaii. By using a 1.148" wide slit and a 2x2 binning, spectra 
with a resolution of R « 29 000 were obtained. These data were 
obtained when the HIRES detector had only a single CCD chip. 
The data are now publicly available and we gathered science and 
associated calibration files from the Keck Observatory Archival 

The data reduction was carried out using the LONG and 
ECHELLE contexts of the MIDAS0 software. Because the bias 



4 Available at http : //koa. ipac . caltech. edu 

3 ESO-MIDAS is the acronym for the European Southern 
Observatory Munich Image Data Analysis System which is developed 
and maintained by the European Southern Observatory. 



level is changing on a time scale of a few minutes, we used the 
over-scan region to compensate for the observed variations and 
to bring all the raw data to an effectively homogenised bias level. 
Master calibration frames were created by averaging the rele- 
vant frames obtained close in time to the science observations. 
The data were then bias, dark and background-illumination sub- 
tracted using standard procedures of the ECHELLE context. The 
orders were traced directly on the science images and a 5-pixel 
window was used to extract the object spectra. Sky spectra were 
extracted from two smaller windows on both sides of the ob- 
ject window. Flat-field and wavelength calibration spectra were 
extracted using the exact same windows as the one used for ob- 
ject and sky extraction. The science and sky spectra were then 
flat-fielded and sky subtracted. Finally, because of little or no 
overlap between the orders, the wavelength calibration was per- 
formed individually for each order. In total, 27 orders were ob- 
served yielding an effective wavelength coverage from 4670 to 
7180 A, although with some gaps between the orders. 



Only two of these four stars could be analysed. The re- 
duced spectrum for MACHO-1997-BLG-045S was not of suffi- 
cient quality to allow for any measurements of equivalent widths 
or line synthesis necessary for a proper abundance analysis, 
and MACHO- 1999-BLG-00 IS appears to be a spectroscopic bi- 
nary. Also, as MACHO- 1999-BLG-006S turned out to be a low- 
luminosity giant after our re-analysis (log g « 2 — 3), the results 
for this star will be presented together with the other similar low- 
luminosity giant stars observed at ESO in a subsequent study. 
The light curve for MACHO-1997-BLG-022S is shown in Fig.Q] 
and event data given in Table Q] 
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Fig. 2. Positions and radial velocities for the 15 
microlensed stars. The arrows represent mea- 
sured radial velocities and one degree cor- 
responds to 70 kms -1 . Upward pointing ar- 
rows indicate positive velocities. The curved 
line shows the outline of the southern Bulge 
based on observation s with the COBE satellite 
jWeilandet al]|1994) . The location of Baade's 
window (BW) is marked by the larger open cir- 
cle. Stars have been labelled with the last num- 
ber in their IDs. 



3. Bulge membership 

3.1. Positions on the sky and radial velocities 

The locations on the sky of the events are shown in Fig. [2] All the 
stars have negative Galactic latitudes because OGLE and MOA 
currently only monitor fields with b < 0. The angular distances 
to the Galactic plane are similar for all events, between 2-5°. 
The measured radial velocities for the stars have been indicated 
in Fig. [2] Arrows upward means positive radial velocities, ar- 
rows downward negative velocities, and the scale in the figure 
is 70 kms -1 per degree. The high variation in v r for the mi- 
crolensed stars is consistent with the high velocity dispersion 
seen for Bulge giant stars (compare, e.g., the rece nt BRAVA ra- 
dial v elocity survey of red giants in the Bulge by Howard et al. 
120081) . 

3.2. Microlensing toward the Galactic bulge 

Because we need to observe the microlensing events wherever 
they occur in the central regions of the Galaxy, we cannot choose 
stars along, e.g., the minor axis to maximise the contribution 
of the Bulge, leading to possible confusion about whether these 
stars are Bulge, disc or halo stars. 

Our current approach to this is to regard the division into 
Bulge and disc for stars within 1 kpc of the Galactic centre as 
a semantic division. There is no evidence for a cold rotating 
extended disc that close to t he Galactic centre for the fields 
that have been studied so far (Ho ward et al.ll2009l) . and the cre- 
ation of the Bulge from the Galactic disc is one of the scenar- 
ios we wish to test. So the question then becomes whether the 
microlensed dwarf and subgiant stars are located in the Bulge 
region, or in the disc on either this s ide or the far side of the 
Bulge. iNair & Miralda-Escudel d 19991) estimate that about 15 % 
of the events toward the Bulge could have source stars belonging 
to the far side of disc, more than 3 kpc away from the Galactic 
centre. On the other hand, more recent theoretical calculations 
of the distance to microlensed sources, assuming a constant disc 
density and an exponential bulge, show that the distance to the 
sources is strongly peaked in the Bulge, with the probability of 
having D < 7 kpc very small dKane & Sahull2006l) . 

Another argument that these are Bulge stars, rather than disc 
stars, are th e large radial veloc ities for stars close to the Galactic 
centre (e.g. lEpstein et al.ll2009l) as well as the fact that there are 



stars with radial velocities in opposite direction s on the same 
side of the Galactic centre (e.g. JCohen et al.|[2008l) . 

In addition, we have prior knowledge of the colours and 
magnitudes of the source stars when unmagnified. The stars we 
observed were all identified as dwarf or subgiant stars at the dis- 
tance of the Bulge based on OGLE or MOA instrumental colours 
and magnitudes, estimated from the offsets from the red clump 
stars in the same field. The parameters for the stars determined 
in this manner have been repeatedly tested against the param- 
eters derived spectrosco pically and the overall co nsistency be- 
tween these results (e.g. lJohnson et al.1 l2007. 2008) again shows 
that these stars are likely located at the distance of the Bulge 
and not the near or far disc. The fact that the latitudes of the 
stars are closer to the Galactic plane may make disc contamina- 
tion more likely for unmagnified sources; however, the combina- 
tion of kinematics, colour-magnitude diagrams and microlensing 
statistics indicate that we are studying a stellar population be- 
longing in the Bulge. 

4. Analysis 

4. 1 . Stellar parameters and elemental abundances 

The determination of stellar parameters and calculation of ele- 
mental abundances w ere carried out as described in method 1 of 
iBensbv et al.l {2009b). Briefly, this method is based on equiva- 
lent width measurements (Wa) and one-dimensional LTE model 
stellar atmospheres calculated with the Uppsala MARCS code 
dGustafsson et alj|1975t lEdvardsson et alj|1993t lAsplund et ail 
119971) . The spectral l i ne list is an expanded version of the list 
used bv lBensbv et al.l(l2003Ll2005l) and is in full given in Bensby 
et al. (in prep.). Equivalent widths were measured using the 
IRAlQ task splot. Gaussian line profiles were fitted to the ob- 
served lines, but in special cases of strong Mg, Ca, Si and Ba 
lines, Voigt profiles were used to better account for the extended 
wing profiles of these lines. 

The effective temperature (T e ff) is found by requiring excita- 
tion balance of abundances from Fe i lines, the microturbulence 
parameter (£ t ) by requiring zero slope in the graph where abun- 
dances from Fei lines are plotted versus the reduced strength 

6 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under co-operative agreement with 
the National Science Foundation. 
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Fig. 3. Diagnostic plots showing absolute Fe abundances versus reduced line strength and lower excitation potential. Open circles 
indicate abundances from Fe i lines and filled circles from Fe n lines. Similar plots for MOA - 2006-B LG-099S, OGLE-2006-BLG- 
265S, OGLE-2007-BLG-349S, and OGLE-2008-BLG-209S can be found in lBensbv etaP d2009bl) . Note the limited number of 
lines for MOA-2009-BLG-259S due to that this star was observed when only the UVES blue CCD was available. 



(log(H^i//l)) of the spectral lines, and the surface gravity (logg) 
from ionisation balance, i.e., requiring that the average abun- 
dances from Fe i and Fe n lines are equal. Only Fe i and Fe n lines 
with measured equivalent widths smaller than 90 mA are used in 
the determination of the stellar parameters. Figure [3] shows the 
diagnostic plots, log e(Fe) versus XogiW^I 'X) and lower excitation 
potential 0f e ), for the stars. 



To relate the elemental abundances to those in the Sun we 
determine our own solar abundances. The equivalent widths we 
measure in the solar spectrum that was obtained by observing 
the asteroid Pallas with UVES on April 11, 2009, show very 
good agreement with the equivalent widths of several solar spec- 
tra (average of Ganymede, Ceres, Vesta, Moon, and sky spectra) 
in Bensby et al. (in prep.). On average the measurements of the 
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Table 2. Stellar parameters, ages, and radial velocities for the sample of microlensed dwarf stars. 1 
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Object 


Teg 


log 9 


ft 


[Fe/H] 


Afoi, Afo.ii 


Age 


u r 




[K] 


[cgs] 


[kms- 1 ] 






[Gyr] 


[kms -1 ] 



OGLE-2009-BLG-076S 


5877 


±96 


4.30 


±0. 


18 


1. 


.61 


±0.16 


-0.72 


±0.07 


66. 


7 


11.7 


± 1 


.9 


+ 128.7 


MOA-2009-BLG-493S 


5457 


±98 


4.50 


±0. 


22 


0. 


.83 


±0.30 


-0.70 


±0.14 


80, 


5 


9.1 


±4 


.0 


-14.5 


MOA-2009-BLG-133S 


5597 


±92 


4.40 


±0. 


27 


1. 


.15 


±0.27 


-0.64 


±0.17 


68, 


7 


9.4 


±4 


.0 


+91.6 


MOA-2009-BLG-475S 


5843 


± 173 


4.40 


±0. 


30 


1. 


.31 


±0.31 


-0.54 


±0.17 


53, 


4 


9.1 


±3 


.7 


+ 137.8 


MACHO-1999-BLG-022S 


5650 


± 106 


4.05 


±0. 


20 


0. 


.30 


±0.45 


-0.49 


±0.14 


97, 


10 


11.7 


± 2 


.3 


+37.6 


OGLE-2008-BLG-209S 


5243 


±65 


3.82 


±0. 


16 


1. 


.01 


±0.11 


-0.30 


±0.06 


146, 


19 


9.5 


±3 


.1 


-173.6 


MOA-2009-BLG-489S 


5634 


±89 


4.30 


±0. 


18 


0. 


.68 


±0.17 


-0.18 


±0.11 


114, 


15 


11.0 


±2 


.3 


+96.5 


MOA-2009-BLG-456S 


5700 


±93 


4.24 


±0. 


18 


1. 


00 


±0.20 


+0.12 


±0.09 


91, 


14 


7.9 


± 1 


.6 


-164.6 


OGLE-2007-BLG-514S 


5644 


± 130 


4.10 


±0. 


28 


1. 


.55 


±0.29 


+0.27 


±0.09 


49. 


8 


6.9 


± 1 


.6 


+ 158.8 


MOA-2009-BLG-259S 


5000 


±400 


4.50 


±0. 


50 


0. 


.50 


±0.75 


+0.33 


±0.40 


16, 


3 


9.4 


±4 


.0 


+83.2 


MOA-2008-BLG-311S 


5944 


±68 


4.40 


±0. 


13 


1. 


.17 


±0.12 


+0.36 


±0.07 


118, 


16 


2.3 


+ 1 


.0 


-34.1 


MOA-2008-BLG-310S 


5704 


±65 


4.30 


±0. 


12 


1. 


.05 


±0.11 


+0.42 


±0.08 


122, 


20 


4.5 


± 1 


.0 


+77.5 


OGLE-2007-BLG-349S 


5229 


±63 


4.18 


±0. 


13 


0, 


.78 


±0.13 


+0.42 


±0.08 


103, 


18 


13.6 


± 1 


.0 


+ 113.0 


MOA-2006-BLG-099S 


5741 


±87 


4.47 


±0. 


15 


0. 


.84 


±0.14 


+0.44 


±0.10 


119, 


21 


2.8 


+ 1 


.3 


+99.0 


OGLE-2006-BLG-265S 


5486 


±70 


4.24 


±0. 


15 


1. 


.17 


±0.12 


+0.47 


±0.06 


92, 


14 


7.9 


+ 1 


.3 


-154.0 



f For the four stars analysed in lBensbv et aD <2009bl) : OGLE-2008-BLG-209S, M OA-2006-BLG-099S , OGLE-2006-BLG-265S, and OGLE-2007-BLG-349S, we 
have here updated the estimations of the uncertainties using the method outlined in Epstein et al. 1 2009). 



Pallas solar spectrum is only 0.3 % larger, which is truly negli- 
gible. Hence, to ensure that the normalised abundances for the 
microlensed dwarf stars are on the same baseline as the sample 
of ~ 700 thin and thick disc dwarf stars in Bensby et al. (in prep.) 
we use the average equivalent widths based on measurements in 
all solar spectra (see Bensby et al. in prep.). 

Final abundances are normalised on a line-by-line basis and 
then we take the median value for each element. In a few cases 
when the equivalent width of an Fe line in the Sun was larger 
than 90 mA, or when a Ti or Cr line were larger than 110 mA, 
and these lines were measured in the Bulge dwarf star, we nor- 
malised the abundance for that line with the average abundance 
from all other lines that were measured in the solar spectrum for 
that element. These cases are marked by "av" in col. 7 in Table|4] 

Final stellar parameters for our targets are given in Table [2] 
All measured equivalent widths and elemental abundances for 
individual spectral lines are given in Table|4] while Table[5]gives 
the normalised abundance ratios. 



4.2. Additional dwarf stars 

In order to increase the sample size of microlensed dwarf 
stars, we include the two st ars MOA-20q9-BLG -310S, and - 
31 IS, recently publ ished bv ICohen et alJ (l2009bl) . and OGLE- 
2007-BLG-514S bv lEpstein et all (120091). The metallicities that 
were found for the se three stars by ICohen et al.l (l2009bl) and 
lEpstein et alJ (120091) are [Fe/H] = +0.41, [Fe/H] = +0.26, and 
[Fe/H] = +0.33, respectively. 

The spectra for these stars were kindly provided by the au- 
thors and we have re-analysed them using our methods in or- 
der to have all 15 microlensed dwarf stars on the same base- 
line. The values we find for these stars are listed in Table [2] and 
they are generally in go od agreement with w hat were found in 
ICohen et al.l (l2009bl) and lEpstein etal.1 (120091) . The main differ- 
ences are that we derive a 240 K higher and 0.2 dex higher 
log g for MOA-2009-BLG-311S, and a 0.2 dex lower logg for 
OGLE-2007-BLG-5 14S. The other differences are within the es- 
timated uncertainties. 



4.3. Error analysis 

A rigorous error analysis as outlined in Eps tein et al.l (120091) has 
been performed for the microlensed dwarf stars. This method 
takes into account the uncertainties in the four observables that 
were used to find the stellar parameters, i.e. the uncertainty of 
the slope in the graph of Fe i abundances versus lower excita- 
tion potential; the uncertainty of the slope in the graph of Fei 
abundances versus line strength; the uncertainty in the difference 
between Fe i and Fe n abundances; and the uncertainty in the dif- 
ference between input and output metallicities. The method also 
accounts for abundance spreads (line-to-line scatter) as well as 
how the abundances for each element reacts to changes in the 
stellar parameters. 

The resulting errors in the stellar parameters are given to- 
gether with the best fit values of the stellar parameters in Table|2] 
The errors in the abundance ratios are given in Table [5] 

4.4. Stellar ages 

Stellar ages were determined as described in Mel endez et al.l 
(2009). We interpo lated a fine grid of q-e nhanced Yonsei-Yale 
(Y 2 ) isochrones by Demarau e et alJ (120041) . adopting [a/Fe] = 
for [Fe/H] > 0, [a/Fe] = -0.3 x [Fe/H] for -1 < [Fe/H] < 0, 
and [a/Fe] — +0.3 for [Fe/H] < -1. At a given metallicity, we 
searched for all solutions allowed by the error bars in T e ff, log g 
and [Fe/H], adopting as final result the median age and as error 
the standard deviation. Column 6 in Table|2]give the median age, 
and column 7 the error. Figure|4]shows the 15 microlensed dwarf 
stars together with the Y 2 isochrones in the log T e ff - log g plane. 

4.5. Check 1: temperatures from microlensing techniques 

De-reddened colours and magnitudes of the sources can be es- 
timat ed using standard microlensing techniques (e.g. lYoo et all 
2004). The method for determining the colour does not make 
any assumption about the absolute reddening, nor about the ra- 
tio of selective to total extinction. It only assumes that the red- 
dening toward the microlensed source is the same as the red- 
dening toward the red clump, and that the red clump in the 
Bulge has (V - I) = 1.05 and 7 = 14.32 re.g.. lJohnson et all 
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Fig. 4. Illustration of the estimation of stellar ages using the a-enhanced isochrones from iDemarque et alj d2004l) . Each set of 
isochrones have been calculated with the same metallicity and a-enhancement as derived for the stars. In each plot the solid lines 
represent isochrones with ages of 5, 10, and 15 Gyr (from left to right). Dotted lines are isochrones in steps of 1 Gyr, ranging from 
0. 1 Gyr to 20 Gyr. Error bars represent the uncertainties in r e ff and log g as given in Table [2] 



Table 4. Measured equivalent widths and calculated elemental abundances for each star. 1 





A 






Sun 


ob09076 


mb09133 


mb09456 


mb()9475 


mb09489 


mb09493 


mb99022 


El. 


[A] 


[eV] 


W.i,e 


e(X) flg 


W A e(X) 


W A e(X) 


W A e(X) 




W A e(X) 


W A e(X) 


W A e(X) 


Ah 
Ah 


5557.06 
6696.02 


3.14 
3.14 


10.4 
44.9 


6.44 
6.62 

















' For each line we give the log gf value, lower excitation potential (xi), measured equivalent widths (W A ), derived absolute abundance (log e(X)). The table is only 
available in the online version of the paper and in electronic form at the CDS via anonymous ftp to cdsarc .u-strasbg . fr (138 . 79 . 125 . 5) or via 
http: //cdsweb. u-strasbg. fr/Abstract .html. 
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Table 3. Comparison of colours and effective temperatures as 
determined from spectroscopy and microlensing techniques ' . 



Object 



M, (V-I)o T 



pi hot 
eff 



(v-iyr 



OGLE 
MOA- 
MOA- 
MOA- 
MACHO- 
OGLE 
MOA- 
MOA- 
OGLE 
MOA- 
MOA- 
MOA- 
OGLE 
MOA- 
OGLE 



2009- 
2009- 
2009- 
2009- 
1999- 
2008- 
2009- 
2009- 
2007- 
2009- 
2008- 
2008- 
2007- 
2006- 
2006- 



BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 
BLG 



076S 
493S 
133S 
475S 
022S 
209S 
489S 
456S 
514S 
259S 
311S 
310S 
349S 
099S 
265S 



4.19 
3.34 
4.19 
4.25 

2.52 
3.39 
2.76 
4.68 
2.91 
3.85 
3.46 
4.54 
3.81 
3.59 



0.67 
0.79 
0.68 
0.59 

0.71 
0.85 
0.66 
0.70 
0.79 
0.66 
0.69 
0.78 
0.74 
0.68 



5775 
5325 
5730 
6150 

5670 
5200 
5870 
5760 
5450 
5880 
5780 
5490 
5620 
5840 



5877 
5457 
5597 
5843 
5650 
5243 
5634 
5700 
5644 
5000 
5944 
5704 
5229 
5741 
5486 



0.65 
0.75 
0.71 
0.65 
0.71 
0.82 
0.71 
0.71 
0.73 
0.97 
0.65 
0.72 
0.87 
0.70 
0.78 



1 Columns 2 and 3 give the absolute dereddened magnitudes and colours based 
on microlensing techniques; based on the colours in col. 3; col. 4 gives the 
inferred effective temperatures using the colour-[Fe/H]-rj.(f calibrations by 
Ramirez & Melendez 1 2005); col. 5 gives our spectroscopic temperatures (same 
as in col. 2 in Table 2); and col. 6 gives the colours, based on the 
colour-iFe/Hl-T^jf calibrations by Ramirez & Melendez 1 2005), that the 
spectroscopic temperatures in col. 5 gives. 



2008; Epstejn_et_alj[2009). The absolute de-reddened magnitude 
and colour are then derived from the offsets between the mi- 
crolensing source and the red clump in the instrumental colour- 
magnitude diagram (CMD). The absolute de-reddened magni- 
tudes and colours for 14 of the 15 microlensed stars are given in 
Tabled Photometry for MACHO- 1999-BLG-022S could not be 
recovered at this time. 

From the c o lour- [Fe/H]-r e ff calibrations by 
Ramirez & Melendez (2005) we check what temperature 



we should expect given the de-reddened colour and the metal- 
licity we determined. On average we find that the spectroscopic 
temperatures are 103 K lower than the ones based on the 
colour-[Fe/H]-r e ff relationships. The top panel of Fig. [5] shows 
a comparison between the two as a function of [Fe/H]. No 
obvious trends can be seen. 

It is also possible to use the the same calibrations by 
Rami rez & Melendez! J2005) to see what (V-I) colours the spec- 
troscopic effective temperatures and metallicities would give. 
These are listed in the last column of Table [3] and the compar- 
ison between photometric and "spectroscopic" (V - I) colours 
are shown in the bottom panel of Fig. [5] On average the spec- 
troscopic colours are 0.03 mag higher, with no discernible trend 
with metallicity. 

The offset that we see between spectroscopic and photo- 
metric values could be a result of the assumed magnitudes and 
colours of the red clump in the Bulge. Previously, it was assumed 
that the red clump stars in the Bulge had the same colour as the 
red clump stars in the Solar neighbourhood ((V - I)o = 1.00). 
Based on the first microlensing events of dwarf stars in the Bulge 
dJohnson et alJl2008t ICohen et al.ll2008h . and additional obser- 
vational evidence dEpstein et a l. 2009), this value was revised to 
(V - /)o = 1.05. Assuming that the spectroscopic temperatures 
are the correct ones, our results indicate that the (V - /)o colour 
of the red clump stars in the Bulge should be revised upwards by 
an additional few hundredths of a dex to (V - I)o = 1.08. 
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Fig. 5. Top panel shows a comparison of our spectroscopic 
T e g:s to those implied from t he col our-[Fe/H]-r e ff calibra- 
tions by [Ramirez & Melendez (2005) versus [Fe/H]. Bottom 
panel shows a comparison between the colours that our spec- 
tro scopic r e n~:s imply (using the colour-[Fe/H]-r e ff calibrations 
by iRamfrez & Melendez! 120051) and the colours based on mi- 
crolensing techniques. Differences are in both cases given as 
spectroscopic — photometric. 
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Fig. 6. Effective temperature and surface gravity versus [Fe/H] . 
In each figure the regression parameters (slope and constant) are 
given, as well as their uncertainties. The error bars in the stel- 
lar parameters represent the total uncertainty (see Table [2] and 
Sect.O- 



4.6. Check 2: trends with metallicity 

Figure [6] shows how T e g and \ogg vary with metallicity. 
Regression lines, taking the errors in both x and y directions 
into account, are shown, as well as the regression parameters and 
their uncertainties. No significant trends with metallicity can be 
seen. 



4. 7. Check 3: signal-to-noise and continuum bias 

The two papers by iBensbv et al.l (2009a b) found the two first 
microlensed dwarf stars with sub-solar iron abundances. In this 
study we find an additional 5 stars with sub-solar [Fe/H]. It has 
been suggested that the reason for this could be that these spec- 
tra have, on average, lower S /N than the other events studied, 
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especially compared to lCohen et alJd20 09b) which all have very 
high [Fe/H]. The lower S /N should then result in that the contin- 
uum was set too low and thus the Wa& underestimated, resulting 
in too low [Fe/H] . 

We looked in greater detail at MOA-2009-BLG-475S, the 
dwarf star with the spectrum that has the lowest S /N, and tested 
if we could make the star more metal-rich. This experiment was 
done by assuming that the star is actually metal-rich and that 
there are many weak lines that makes it difficult to identify the 
level of the continuum. The best way to set the continuum is then 
to assume that the high points in the spectrum are indeed the 
continuum. We re-measured the star under this assumption and 
determined new stellar parameters. On average the equivalent 
widths became 10 mA larger. This resulted in a change of the 
metallicity of +0. 1 dex, but the other stellar parameters (log g, 
r e ff, and £ t ) were intact and did not change. Hence, we find it 
unlikely that our low-metallicity stars could be high-metallicity 
stars resulting from an erroneous analysis of their relatively low 
S IN spectra. 

4.8. In summary 

Our analysis shows that the eight microlensed stars that we ob- 
served were successfully selected to be dwarf stars, varying in 
metallicity from [Fe/H] = -0.72 to +0.5 4. The first results for 
OGLE-2009-BLG-076S were presented in lBensbv etall ((2009a) 
where we found a metallicity of [Fe/H] = -0.76. The refined 
analysis in this study gives a metallicity of [Fe/H] = -0.72, and 
it still holds (barely) the place as the most metal-poor dwarf star 
in the Bulge, with MOA-2009-BLG-133S just 0.01 dex higher. 

The uncertainties in the stellar parameters are below or 
around 100 K in T e ff, around 0.2 dex in \ogg, and around 0.1- 
0.2 dex in [Fe/H] (see Table |2j. The clear exception is MOA- 
2009-BLG-259S where errors are exceptionally large. This is 
due to difficulties arising from the very limited wavelength cov- 
erage of the UVES spectrum that was obtained when only the 
blue CCD was available. The fact that it also turned out to be a 
metal-rich star at [Fe/H] = +0.32, further increases the errors 
due to line blending and uncertainties in the placement of the 
continuum. Also, the lack of weak lines, due to the high [Fe/H], 
made it especially difficult to determine the microturbulence pa- 
rameter (see Fig. [3). A spectrum covering the whole optical re- 
gion of MOA-2009-BLG-259S was obtained by another group 
using the HIRES spectrogra ph, and they find a 0.2 dex higher 
metallicity than what we do dCohen et alj |2009a). Although the 
metallicity is in reasonable agreement with what others have 
found we think that the errors in T e g and log g are so large that 
this star is not conveying any information in the [Fe/H], age, or 
abundance plots. Therefore we do not include this star in the fol- 
lowing discussions. 



5. Metallicity distributions of dwarfs and giants 

The most recent spectroscopic study of a large homo g eneou s 
sample of giant stars in the Bulge is by IZoccali et al.l (|2008). 
Using FLAMES, the multi-fibre spectrograph at the Very Large 
Telescope, they studied a sample of 521 giant stars at three lat- 
itudes in the Bulge: 204 stars in Baade's window at I ss -4°; 
213 stars at I * -6°; and 104 stars at I * -12°. The 15 mi- 
crolensed dwarf stars observed so far are all located at similar 
angular distances from the Galactic centre as Baade's window 
(see Fig. H}. Therefore, o nly the 204 stars in the / = -4° field 
from[Zoccali et al. (2008) will be used for comparison. 
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Fig. 7. The top panel shows the MDF for the 204 giant stars in 
Baade's window from lZoccali et al.l(l2008l) . and the middle panel 
the MDF for the fourteen microlensed dwarf stars (MOA-2009- 
BLG-259S excluded). Their cumulative metallicity distributions 
(CMDF) are shown in the bottom panel (giant CMDF marked by 
solid line, and the dwarf star CMDF by dash-dotted line). The 
two-sample Kolmogorov-Smirnov D statistic (maximum verti- 
cal distance between two distributions) and the corresponding 
significance level, prob, of D are indicated. 



The average metallicity of the 14 microlensed dwarf and sub- 
giant stars in the Bulge (MOA-2009-BLG-259S excluded, see 
Sect.gS) is ([Fe/H]) = -0.08 + 0.47. This is in agreement with 
the average metallicity of the 204 RGB stars in Baade's window 
that have ([Fe/H]) = -0.04 + 0.40. However, when comparing 
the two distributions, a two-sided Kolmogorov-Smirnow (KS) 
test gives a significance level of the null hypothesis, that they are 
drawn from the same distribution, of 30 % (see Fig. |7J. Hence, 
we can not reject the null hypothesis that the MDF for the 14 
microlens ed dwarf stars and the MDF for the 204 Bulge RGB 
stars from lZoccali et al. I J2008h are identical. 

As discussed by, e.g., ISantos et al.l (120091) it is possible that 
a systematic shift in [Fe/H] between analyses of dwarf stars and 
giant stars, perhaps by as much as 0.2 dex, exists. However, the 
difference in the average metallicity between dwarf and giant 
stars is essentially zero. If there was a real difference of 0.2 dex 
between the dwarf stars and the giant stars, how many dwarf 
stars would we need to observe in order to statistically measure 
that difference? 

In order to estimate the number of stars required to reject 
the null hypothesis, that the MDF for the sample of 204 RGB 
stars, S rgb , is not different from the MDF for the observed dwarf 
stars, we will use a non-parametric bootstrap method. First, we 
construct a sample, Srgb.o.i, of stars that is shifted 0.2 dex from 
the original RGB sample, S rgb- Secondly, we bootstrap n num- 
ber of stars from sample Srgb,o.2, mus creating a new sample, 
Sn,o.2, that contains n stars. Next, we perform a two-sided KS- 
test between samples S rgb and S n ,o.2- If me KS-test yields that 
the distributions are not the same at the 95 % confidence level, 
we reject the null hypothesis. We repeat this process i - 10000 
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Fig. 8. Upper panel: Probability of detecting a difference be- 
tween the Bulge RGB MDF and the microlensed dwarf MDF 
as a function of the number of stars in the dwarf sample. Solid 
and dashed line indicate the probability if the difference is 0.2 
and 0.1 dex, respectively, in the mean of the MDFs. Dotted line 
indicates the probability if the difference is that the microlensed 
dwarf MDF has more metal rich, [Fe/H] > 0.0, stars (i.e, a more 
symmetric MDF) than the RGB MDF. Lower panel: The signif- 
icance of the agreement between the average of the metallicities 
of the RGB MDF and microlensed dwarf MDF, as a function of 
a speculated (still undetected) difference. 



times and take the average, (p), of p l7 where p\ — 1 if the null hy- 
pothesis was rejected, and otherwise (type II error). Thus, (p) 
is our probability to identify an intrinsic difference of 0.2 dex 
in the mean in the MDF for RGB and microlensed dwarf stars. 
Figure [8] shows how this probability varies with the number of 
stars in the hypothetical dwarf star sample. In order to statisti- 
cally verify a difference of at least 0.2 dex, at the 95 % level, we 
need to observe around 40 stars. Also in Fig. [8] we show how 
the probability varies if we want to verify a difference of only 
0.1 dex between dwarfs and giants. Note that the detection of 
such a small difference would require the observation of many 
more microlensed dwarf stars. 

Above we shifted the entire MDF by 0. 1 and 0.2 dex and did 
not consider other statistical parameters that describe the MDF 
such as variance, skewness and kurtosis. However, these param- 
eters will most likely make it easier to reject the null hypothesis, 
that the distributions are the same, if they are considered in the 
test. Additionally, we only considered 0. 1 and 0.2 dex as a differ- 
ence between the samples. A larger difference will also make it 
easier to reject the null hypothesis. Our estimate of (p) is there- 
fore a lower limit. 

Additionally, the average metallicity for the 14 microlensed 
dwarf stars is in agreement with the average metallicity of the 
RGB stars. The question is then, what is the significance of this 
result, given the low number statistics of the microlensed dwarf 
stars, i.e., can we rule out a shift of 0.2 dex between the MDF 
for the RGB stars and the MDF for the dwarf stars? We estimate 



the significance of the agreement between the average metal- 
licities using 10 5 Monte Carlo realisations. In each realisation, 
we draw 14 stars randomly from the RGB sample. We deter- 
mine the significance of the agreement between the average of 
the metallicities by computing the fraction of realisations that 
fail to produce a dwarf star average metallicity lower or equal 
to the average metallicity of the RGB stars. Next, we shift the 
MDF for the RGB stars by a small positive amount, 6, and re- 
peat the above given exercise. Figure [8] lower panel, shows the 
significance as a function of 6. We find that for 5 ~ 0.2, the sig- 
nificance has dropped to 0.05 which indicates that we can rule 
out the possibility that there is a systematic shift in [Fe/H] of 
0.2 dex between the MDF for the RGB stars and the MDF for 
the dwarf stars. However, we need to put this in relation to the 
standard deviation of the difference between the averages given 
by 

o- d = -vp + (1) 

where subscript 1 and 2 indicate the RGB star sample and dwarf 
star sample, respectively, cr is the sample standard deviation, and 
n is the number of elements in each sample. We find a <x rf of 0. 12, 
which is not even 2cr from 0.2 dex. Thus, given the small num- 
ber of dwarf stars in our sample, «2 = 14, we can not rule out 
a possible systematic shift of 0.2 dex between the MDF for the 
RGB stars and the MDF for the dwarf stars, even though the dis- 
tributions at this point have the same average metallicity. Thus, 
more microlensed events are needed for the agreement between 
the averages to be significant. 

The above discussion focused on a systematic offset between 
the dwarf and giant abundances, caused either by systematic ef- 
fects in the data analysis or differences in the surface composi- 
tions of dwarfs and giants because of their differing evolution- 
ary states. However, instead of an overall shift, the dwarf MDF 
may be different from the giant MDF only at the high metallic- 
ity end. Kalira Tet al.1 (120071) suggested that up to 40 per cent of 
the stars at the metallicity of NGC6791 ([Fe/H] = +0.3) skip 
the He burning phase, resulting in a depletion of the HB and 
AGB phases. Therefore, an MDF based on giants may not re- 
flect the MDF for the dwarf stars (i.e. there are metal rich stars 
missing in the MDF for the RGB stars). The question is then, 
how many microlensed dwarf stars are required to reject the null 
hypothesis that the MDF for these dwarfs, that has an excess 
of metal-rich stars compared to Srgb, is no different than the 
MDF for S rgb? Based on this, we construct a dwarf star sam- 
ple with more metal-rich stars than Srgb, 5 , RGB,Kaiirai, under the 
assumption that 100% of the dwarfs with [Fe/H]< 0.0 evolve 
to RGBs, but that the number of dwarfs that evolve to RGBs 
decrease linearly down to 60% at [Fe/H] = +0.3. We extrap- 
olate this linearly for more metal-rich stars, i.e. going down to 
~34% at [Fe/H] = +0.5. Thus, we add stars, more metal-rich 
than [Fe/H] = 0.0, to S rgb to create a hypothetical sample of 
dwarf stars in the Bulge, SRGB.Kaiirai- Additionally, we add stars, 
drawn from a Gaussian distribution centred on [Fe/H] = +0.2 
with cr = 0.4, to S RGB.Kaiirai for the metal rich region where there 
are no observed RGB stars. Thus, making the sample more sym- 
metric. The dotted line in Fig. [8] shows how (p) varies for this 
analysis with the number of stars in SRGB.Kaiirai- We note that 
about 100 stars are required in order to verify that S rgb is dif- 
ferent than S RGB.Kaiirai- 

Also, there has been some claims that the microlensing event 
itself alters the spectrum of the source star, and that this is the 
reason for the very high metallicities of some of the dwarf stars 
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Fig. 9. Ages versus [Fe/H] for 14 microlensed dwarf and sub- 
giant stars in the Bulge (MOA-2009-BLG-259S has been ex- 
cluded). The error bars represent the total uncertainty [Fe/H] and 
age (see Tableland Sects. 1431 andl4~4l 



in the Bulge (e.g. IZoccali et alJl20 08). Recently. iJohnson et alJ 
(2009) investigated the effect that differential limb darkening 
has on abundance analysis of microlensed dwarf stars. They do 
find changes in the measured equivalent widths as a result of the 
differential limb darkening. However, the effect is very small, 
leading to changes in T e ff less than 45 K, log g less than 0. 1 dex, 
and [Fe/H] less than 0.03 dex. Hence, a possible differential limb 
darkening can not be responsible for the MDF discrepancy (if 
any) between dwarf and giant stars. 

In summary, it i s evident that the ex tremely super-metal-rich 
MDF proposed by ICohen et alJ (l2009bl) . exclusively based on 
dwarf stars with super-solar [Fe/H], has shifted toward lower 
metallicities. The MDF of the 14 microlensed dwarf stars is still 
poorly determined, currently being double peaked with excesses 
of low- and high-metallicity stars. Whether this is an effect 
of small number statistics or not is unclear. More microlensed 
events will certainly help to clarify the dwarf star MDF and to 
refine the comparison with the giant star MDF . Also, an o ut- 
standing issue is the puzzle presented in lCohen"et al. (2009a|) of 
the correlation between A max and [Fe/H], which we hope to di- 
agnose as future events are observed. 



6. Ages and metallicities 

The 14 microlensed dwarf and subgiant stars (MOA-2009-BLG- 
259S excluded, see Sect l4.8l ) in the Bulge have an average age of 
8.4+3.3 Gyr. Figure|9]shows the age-metallicity diagram and it is 
evident that stars with sub-solar [Fe/H] all have high ages, while 
at super-solar [Fe/H] there is a large spread in ag e, covering the 
whole age-range see n in the Galactic disc(s) (e.g.. lTwarogl l980: 
Feltzing et al. 2001). That the large age-range seen for super- 
solar [Fe/H] is real is exemplified by two stars: OGLE-2007- 
BLG-349S and MOA-2008-BLG-310S. The first star has a high 
age while the second has a low age. That the ages are robust can 
be seen in Fig. [4] and we also note that their stellar parameters 



have small errors and are derived from a large number of Fe i 
and Fen lines (see Tableland Fig. [3). 

In spite of the small sample, it is notable that we only see the 
low ages for the metal-rich stars while all stars with sub-solar 
[Fe/H] are consistent with the classical view of the Bulge as an 
old population (e.g.. iHoltzman et al.|[l993t iFeltzing & Gilmorel 
2000: IZoccaliet all l2003h . As the stars with sub-solar [Fe/H] 
also have enhanced levels of a-elements (see Fig.fTUb. these stars 
all appear to adhere to the classical picture of the Bulge as a 
stellar population that formed rapidly early in the history of the 
Galaxy (see, e.g., models and discussions in lMatteuccill200 1|) . 

Overall, the evidence for young stars in the Bulge is scarce. 
For instance, the extremely de ep CMD of ~ 180 000 field stars in 
the Bulge bv lSahu e t al. (2006) show no traces of a young popu- 
lation. It is therefore surprising to find three (MOA-2008-BLG- 
310S, MOA-2008-BLG-311S, and MOA-2006-BLG-099S) out 
of 14 stars to have young ages. At this point we can only spec- 
ulate on their origin. One interpretation would be that the older 
stars are all bona fide Bulge stars while the young, metal-rich 
stars are disc interlopers. In the G alactic disc a young age and 
a high metallicity is common (e.g.. lTwarog||1980T IFelt zing et al. 
2001; Nor dstrom et al.l l2004'). Also, the innermost Galactic disc 
is ex pected to be more m etal-rich than the Solar neighbourhood 
(e.g. lColavitti et al.i r2009). However, we still do not now if it is 



supposed to be young too. 

It should furthermore be noted that these young Bulge stars 
are not brighter than the main old turnoff, they are just too blue to 
fall on old isochrones (see Fig. |4j». Also, there are some theoret- 
ical limitations of the isochrone fitting method. First, isochrones 
at [Fe/H] > +0.3 have very few calibrators, and, second, the 
colour of the main sequence is strongly affected by the Y (he- 
lium) content, which for the Bulge, or any population with such 
high metallicity, is poorly known. However, we find a whole 
range of ages at high metallicities, so we don't see a bias in our 
ages. 



7. Abundance trends in the Bulge 

7. 1 . General appearance 

Figure [10] shows the abundance results for 13 of the 15 the 
microlensed dwarf and subgiant stars in the Bulge (MOA- 
2009-BLG-259S w ere excluded, see Sect. EU and for the 
lEpstein et all (2009) star OGLE-2007-BLG-514S we only de- 
termined stellar parameters and Fe abundances.). 

Regarding the <?-elements (Mg, Si, Ca, and Ti), the Bulge 
dwarfs show enhanced [a/Fe] ratios at sub-solar [Fe/H], that de- 
cline when approaching solar metallicities. At higher metallic- 
ities the [a/Fe] ratios are around or slightly higher than solar. 
The oxygen trend is similar to the a-element trends at sub-solar 
metallicities, but differs at super-solar [Fe/H] where it continues 
to decrease. The oxygen abundances that are based on the in- 
frared triplet lines at 773 nm have been NLTE corrected the a c- 
cording to the empirical formula given in Bensb y et al.l ((2004). 

Generally, abundance trends of the dwarf stars in the Bulge 
are very well-defined. The scatter in the [Ti/Fe]-[Fe/H] plot for 
instance is remarkably low. The canonical interpretation of the 
plateau of high a-element abundances relative to iron at low 
metallicities is due to early and rapid chemical enrichment of 
the Bulge by massive stars. When these stars die, they explode 
as core-collapse supernovae, producing a lot of <?-elements rel- 
ative to iron. At some point low-mass stars start to contribute 
to the chemical enrichment, and since these produce less of the 
of-elements the [a/Fe] ratio will start to decline. 
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Fig. 10. Abundance results for 13 microlensed dwarf and subgiant stars in the Bulge (marked by filled bigger circles). Thick disc 
stars from Bensby et al. (2010, in prep.) are shown as small circles, and the solid line is the running median of the shown thick disc 
sample, and the dashed line the running median of the (not shown) thin disc sample from Bensby et al. (2010, in prep.). The error 
bars represent the total uncertainty in the abundance ratios (see Sect. l4.3l and Table[5]). 
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Fig. 1 1. Comparisons of [Na/Fe], [Mg/Fe], and JAl/Fel between our microlensed dwarf stars (filled circles), giants from 
Fulbright et al.l d2007l) (empty squares), giants from lLecureur et al] d2007l) (asterisks), nearby thick disc stars from Bensby et al. (in 
prep.) (small open circles). Solid and dashed lines represent the running median of the thick and thin disc stars, respectively (same 
as in Fig. [Toll. 



The enhanced [Na/Fe] « +0.1 for metal-r ic h dis c 
stars was already noticed bv lEdvardsson et al.l (1 19931) . 
IFeltzing & Gustafssonl d!998l) . IShi et all d2004h . and also 
in iBensbv et alJ d2005h an upturn in [Na/Fe] can be seen. Our 
microlensed dwarf stars are in full agreement with these disc 
results. 

The [Ba/Fe] trend is flat and slightly under-abundant com- 
pared to the Sun for all [Fe/H]. 



7.2. Comparisons to the Galactic thick disc 

In Fig.[10]we also show the thick disc abundance tr ends based on 
dwarf stars i n the Solar neighbourhood (taken from Bensbv et al. 
120031 120051 and Bensby et al., in prep.). These thick disc stars 
have been analysed using the exact same methods (spectral line 
lists, atomic data, model stellar atmospheres, etc.) that we use 
for the microlensed dwarf stars. Hence, any differences between 
the stars from the two stellar populations should be real, and not 
due to unknown systematic effects. 

The first thing that can be taken away from Fig.[10]is that, at 
sub-solar metallicities, the abundance trends for the Bulge dwarf 
stars are very similar to those of the thick disc stars. The solid 
lines shown in the figures indicate the median abundance ratio 
as a function of metallicity for thick disc stars. Even though the 
appearance of the abundance trends at sub-solar [Fe/H] are very 
similar between the Bulge and the thick disc, it is also evident 
that for many elements the Bulge stars appear to be slightly more 
enhanced than the median thick disc. The Bulge stars seem to oc- 
cupy the upper envelope of the thick disc abundance trends. Part 
of this apparent shift in the abundance trends between the two 
populations could be due to that the thick disc sample is kinemat- 
ically selected, and hence will unavoidably be m ixed, to some 
degre e, with kinematically hot thin disc stars (see Bensbv et al. 
2007). The median thick disc line that we show in Fig. [TOlwill 
then be slightly too low. However, this effect should only be im- 
portant when approaching solar metallicities, and we do see a 
shift between the Bulge and the thick disc at lower [Fe/H] as 
well, where the kinematic confusion between the thin and thick 
discs should be negligible. However, the shift is not for all el- 
ements, and it is small, on the order of 0.05 dex or less. More 
microlensing events will help us to clarify if this shift is real or 
not. 



A possible link between the Bulge and the Galactic thick 
dis c based on simi l arities of abundances was first suggested 
by iProchaska et al.l d2000l) . pointing out the "excellent agree- 
ment" between the abundance ratios in their sample of ten thick 
disc s tars to those of the Bulge giants of iMc William & Rich! 
d!994l) . As the thick disc sample of IProchaska et all d2000l) only 
reached [Fe/H] « -0.4, the downturn in [g/Fe l that we now 
see in the thick disc at [Fe/H] « -0.35 (e.g. IFeltzing et al.l 
120031: IBensbv et al.ll2007h was at that time not known. Hence, 
when later studies of giant stars in the Bulge showed that 
the ra/Fe] remained high even at super-solar metallicities (e.g. 
Fulbright end] [2007b . in contrast to the declining thick disc 
trends, the possible connection between the Bulge and the 
thick disc became le ss clear. It should also be noted that 
Prochask a et al.l d2000l) did n ot analyse both t he Bulge and thick 
disc samples. Also, recently. iMelendez et alJ d200 8) presented a 
consistent analysis of giant stars in both the Bulge and the thick 
disc that foun d a similarity be tween them for C, N, O (recently 
confirmed bv iRvde et alJl2009l) . The agreem ent is extended to 
other g-elements in the upcoming study by lAlves-Brito et al.l 
J20T0I) . However, the first confirmation that the Bulge and the 
thick disc have similar abundanc e patterns based on dwarf stars 
came from Bensbv et al.l d2009blal) . and with this study it now 
appears well established that the Bulge and the thick disc have 
had, to some degree, similar histories. 

Before interpreting the apparent agreement between the 
Bulge and thick disc abundance trends, it is important to recog- 
nise that the MDFs for the thick disc and the Bulge differ in the 
average as well as width of their MDFs. The metall icity distribu- 
tion f or the thick disc peaks at [Fe/H] » -0.6 (e.g JCarollo et alJ 
120091). and its met al-rich tail likely reach solar metallicities 
dBensbv et all2007l) . Stars with [Fe/H] > that can be kinemat- 
ically classified as thick disc stars, are heav ily biased to belon g 
to the high-velocity tail of the thin disc (Bens bv et al.l [2007). 
Even though the average metallicity of the Bulge and the shape 
of its MDF is under debate (see Sect. [5]) it is clear that it spans 
the full range of metallicities of the thick disc, and in addition 
reaches very high, super-solar [Fe/H]. Hence, the agreement, as 
regards the abundance trends, between the Bulge and the thick 
disc can currently only be established for sub-solar metallicities. 
The comparison that can be made at super-solar [Fe/H] is be- 
tween the Bulge stars and thin disc stars (dashed lines in Fig.lTOb 
that reach similar high metallicities. As the abundance trends in 
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Fig- HO] show, the agreement between the Bulge and the thin disc 
at [Fe/H] > is good. Looking at, e.g., the [Na/Fe], [Ni/Fe], and 
[Ba/Fe] trends, which are not simply flat at the highest [Fe/H], 
the Bulge stars nicely extend th e trends that a re see n for the 
nearby thin disc stars. However, lEpstein et alj d2009l) found a 
[Ba/Fe] value noticeably below the disc for OGLE-2009-BLG- 
514S. The Ban lines are strong and very sensitive to the micro- 
turbulence parameter, and therefore measuring the abundances 
of another heavy element, such as La, would be a useful check. 

With this study where we compare Bulge dwarf stars with 
disc d warf stars, and with the recent studies bv lMelendez et al] 
(2008li) and lAlves-Brito et al.l d2010l) . all using internally con- 
sistent methods, it appears clear that the Bulge and thick disc 
abundance trends are similar. The agreement between the Bulge 
and the thick disc means that conclusions from recent theoreti- 
cal works, developed under the assumptions that the abundance 
trends in the Bulge and thick disc are different, may not be valid. 

We caution that the comparisons that we do and that have 
been done by others are between the Bulge stars close to the 
Galactic centre and thin and thick disc stars in the vicinity of the 
Sun, i.e., at a distance of approximately 8kpc from the Galactic 
centre. In order verify a possible connection between the Bulge 
and the Galactic thick (and thin) disc(s) we need disc stellar sam- 
ples much closer to the Bulge, say at 4kpc from the Galactic 
centre. At that distance the contamination of Bulge stars in a 
disc stellar sample should be small, allowing us to directly com- 
pare the populations. If, for instance, the star formation has been 
more rapid in the inner disc than in the Solar neighbourhood, or 
in the outer disc, we should expect the "knee" in the c-element 
abundance trends at higher [Fe/H] than what we see in the Solar 
neighbourhood. Then that could explain the apparent slight shift 
in the abundance trends between the Bulge and the thick disc 
in the Solar neighbourhood that we see. However, no such inner 
disc sample is currently available. 

7.3. Comparisons to Bulge giants 

It should be noted that some recent studies of Bulge giant stars 
claim that the Bulge is enhanced in the a-eleme nts with respect 
to the stars of the Galactic disc. For instance, iRich & O riglia 
compare their Bulge giants to nearby thi n disc giants 
and fi n d the Bulge on e s to b e m ore enhanced ; Zoccali et alJ 
(120061) : IFulbright et all (120071) and iLecureur et all (120071) find 
their Bulge giants to be more enhanced than comparison sam- 
ples of nearby thin and thick disc dw arf stars. In contrast, the 
re cent study bv | Melendez et alJ d2008), and the upcoming study 
by lAlves-Brito et all d2010l) . compare Bulge giants with nearby 
thick disc giants, and find the Bulge and thick disc abundance 
trends to be similar. 

In Fig. Q~T] we compare [Na/Fe], [Mg/Fe] and [Al/Fe] 
between our microl ense d Bulge dwarfs, the giants from 
iFulbright etail d2007l) and ILecureur et all d2007l) . and the Solar 
neighbourhood thin and thick disc dwarfs from Bensby et 
al. (2003, 2005, and in prep.). It is evident that the giants from 
ILecureur et al ] d2007l) are much more enhanced in Na, Mg, and 
Al than any of the other compari son samples, w i th our dwarf 
stars having lower abundances, the Fulbrigh t et al.l d2007|) giants 
having intermediate abundances, and the ILecureur et al.l d2007l) 
giants having the highest abundances. The higher abundances of 
the giant stars might be due to the way these abundances have 
been normalised to the Sun. For dwarf stars with spectral types 
similar to the Sun, the normalisation is straightforward, it is just 
to analyse the Sun in the same way and then subtract the so- 
lar abundance. However, for giants, having very different log g:s 



and r e ff:s than the Sun, the normalisation is usually done to an- 
other standard star, such as, e.g., fi Leo. This could partly explain 
the le vels of the abun d ances , progressing from our dwarf stars, 
to the IFulbright et al.l ([2007) giants, and to the ILecureur et alJ 
d2007l) giants. 

Looking closer at the ILecureur et al .1 d2007l) data it appears 
as if there is a very well defined lower envelope just above the 
thick disc trend, and on top of that a large scatter of stars with 
higher abundance ratios, spreading upwards. We suspect that this 
might be due to line blending, and possibly lack of accurate con- 
tinuum points, in the extremely crowded spectra that metal-rich 
giants have. The lower envelope that is seen in the lLecureur et al.l 
(2007) is what could be expected if blending is the case. A num- 
ber of stars that have less (or very few) blends forming the well- 
behaved lower envelope of the trends, while others are more 
affected, leading to randomly increased equivalent widths, and 
hence randomly increased abundances. 

7.4. No Na-0 anti-correlation 

The Na-O anti-correlation that has seen in all gl obular clusters 
studie d so far is not present for field stars (e.g., ICarretta et al.l 
2009) . Therefore, o ne of the most striking results from 
ILecureur et al.l (2007) was the Na-O anti-correlation that they 
found for their Bulge giants, and they claim it is probably an 
effect of the Bulge chemical evolution. However, Fig. [12] (left 
panel) shows the Na-O plot for our microlensed dwarf stars and 
there is no Na-O anticorrelation present. Instead [Na/Fe] versus 
[O/Fe] is flat. Given the large spreads that lLecureur et all d2007l) 
see for Na, Mg, and Al (see Fig. fTTb it appears likely that their 
results are affected by systematic errors (blending presumably). 



7.5. Metallicity dependent oxygen yields 

In Fig. [12] (right panel) the evolution of [O/Mg] vs [Mg/H] is 
shown, which reveals a very tight correlation. Our results from 
microlensed dwarfs have significantly less scatter than the corre- 
sponding values for Bulge red giants (e.g. Lecureur et al.l l2007: 



Fulbright 'et"al.l2007l: lAi ves-Brito et al.ll2010l) presumably due to 
the lesser influence of blends and better determined stellar pa- 
rameters. The dwarf-based slope is also slightly shallower than 
for the giants. A declining trend in [O/Mg] towards higher metal- 
licity is expected with traditi onal metallicity-dependen t super- 
novae yields such as those of IWooslev & Weaver! d 19951) but our 
slope is steeper, which may signal an metallicity-depen dence in- 
troduced for example by mass- loss in massive stars (e.g.lMaedet 
ll992t | Mevnet & Maederl2002l) as discussed by Mc William et al. 
(2008J) and lCescutti et al.l (I2009h . These findings are based on the 
increased Wolf-Rayet stellar wind efficiency at higher metallic- 
ity that removes a larger fraction of He and C before they are 
converted to O and thus decreasing the O production but leaving 
the Mg-yield largely unaltered. The similarity between the Bulge 
and disc results implies that these metallicity-dependent nucle- 
osynthesis yields are a general feature but also argues against 
substantial differences in the initial mass function (IMF) be- 
tween the two populations. 



8. The origin of the Galactic bulge 

To summarise our observations, we find that the stars with high 
[ff/Fe] and low-metallicity in the Bulge are old, but the stars with 
high metallicities and solar (or subsolar) [a/Fe] ratios span a 
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Fig. 12. [Na/Fe l versus ["O/Fe] (Left p anel) and [O/Mg] versus [Mg/H] (Right panel) for our m icrolensed dwarf stars (filled circles), 
the giants from Fulbrig ht et al.1 (|2007 | ) (open squares), the giants from Lecur eur et alJ |2007) (asterisks), and thick disc dwarf stars 
in the Solar neig hbourhood dBensbv et alJl2003Ll2005L and in prep.) (small open circles). 



range of ages from 10 Gyr to 3 Gyr. What does this mean for 
the origin of the Bulge? 

The old age, high [a/Fe] population can be explained 
through the standard chemical pattern in systems dominated by 
Type II SNe that is common to all stellar systems. That the 
high [a/Fe] ratio persist to [Fe/H] > -0.6 indicates that star 
formation proceeded very efficiently in the event(s) that cre- 
ated the Bulge. Such events could be early mergers of subha- 
los which drive efficient star formation as well as contributin g 
their own high [a/Fe] and low [Fe/H] stars (iRahimi et alJ l2009). 
or early fragmentation of the disc into clu mps of stars and ga s 
which then rapidly merge to form the Bulge (llmmeli et al. 2004). 
These results may also be consistent with the secular evolu- 
tion of the disc, depending on the age-metal licity-fa/Fe] relation 
presen t in the inner disc. In the models of ISchonrich & Binnevl 
(2009), the inner disc is composed of stars that have old ages 
and higher [a/Fe] at the same [Fe/H] than stars formed in the 
thin disc in the Solar neighbourhood. Some of these inner disc 
stars then migrate outwards to form the local thick disc. Thus 
it would not be surprising that the Bulge (=puffed up inner 
disc) stars should be chemically similar to the local thick disc 
(=migrated in ner disc) stars, which is e xactly what we see. 
The model by Scho nrich & Binnevl d2009l) requires a number of 
assumptions/approximations/parameterisation of the migration 
and heating processes, and is tuned to explain the thick disc in 
the Solar neighbourhood, so observational evidence (still miss- 
ing) of the nature of the inner disc is, once agai n, important. 
Our re sults are not in agreement with the models of llmmeh et al.l 
(2004) where the gas cools less efficiently and the instability that 
forms the bar sets in at later times, because those models predict 
a decline in [a/Fe] starting at much lower metallicities (see their 
Figure 10). 

The solar [a/Fe] and high metallicity stars that span a range 
of ages are more of a puzzle. The old stars imply that part of the 
Bulge got a head start on its chemical evolution, so that some 10 
Gyr old stars were formed from a population that was already 



producing Type la SNe, while other 10 Gyr stars were forming 
out of gas that had just been enriched with Type II SNe. Old, 
low [a/Fe] stars are found in the dwarf spheroid galaxies, but 
accreting a Sagittarius-like object would not explain the pattern, 
because th e low [a/Fe] stars in Sagittarius have too low metallic- 
ities (e .g., IVenn et aT1l2004h . In the simulations of Rahimi ^t al.l 
(2009), the bulges experience a series of mergers over a period 
of ~5 Gyr, leading to populations of "old", "intermediate" and 
"young" stars. The distributions of [Mg/Fe] for these populations 
do show some old stars with low [Mg/Fe], (as well the expected 
shift to low [Mg/Fe] for the younger stars) but they caution that 
their code suppresses mixing among gas particles, leading to ar- 
tificially high abundance ratio dispersions. Nonetheless, this is 
what we see in our data, so perhaps the Galaxy found a way 
to suppress mixing as well, maybe with a merg er history dif- 
ferent than the two cases in lRahimi et al. (2009), where two of 
the subclumps that merged had different starting times for star 
formation relative to today. 

Finally, the younger, low [a/Fe] and high metallicity stars 
show that star formation persisted in the components that cre- 
ated the B ulge. These s t ars a re seen in the disc fragmentation 
models of llmmeh et al.l (120041) . where, for the cold gas model 
that agrees with the turnover in the [a/Fe] vs. [Fe/H] diagram, 
30 % of the baryons are not rapidly converted to stars. Instead, 
low level star formation occurs for several Gyr and produces a 
peak in the [Mg/Fe] histogram of the Bulge at [Mg/ Fe] = -0.2. 
This k ind of stars are also seen in the models of Rahimi et al. 
(120091) where new stars formed after mergers at later times are 
polluted with Type la ejecta as well. Finally, depending on the 
star formation of the thick disc, these stars may be present in that 
component and then used to make the Bulge. 

To distinguish further among these models, the inner disc 
of the Galaxy needs to be better characterised observationally. 
In addition, more Bulge dwarfs with accurate ages, metallicities 
and abundance ratios would help clarify whether the age spread 
is confined to the higher metallicities and whether the oldest low 
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[oTFe] stars are as old as the lower metallicity high [a/Fe] stars. 
Finally, the elements produced in Type la SNe are not the only 
chemical evolution "clock" available, and measuring elements 
produced in AGB stars (C, N and s-process), for example, would 
test whether chemical evolution really began earlier for some 
stars now in the Bulge than for others. 

9. Summary 

With this study we have doubled the number statistics on the 
data for microlensed dwarf and subgiant stars in the Bulge. All 
stars have been observed with high-resolution spectrographs and 
from a detailed elemental abundance analysis we present results 
for O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, Y, and Ba. The 
method we utilise is identical to the method used for a large 
sample of 702 F and G dwarf stars in the thin and thick discs 
in the Solar neighbourhood. Therefore, any differences between 
the Bulge stars and the disc stars should be genuine, and not due 
to unknown (systematic) uncertainties. We have also determined 
stellar ages for the stars, a task that is impossible to do with giant 
stars. 

The main results and conclusions that can be drawn from our 
sample of 15 microlensed stars in the Bulge are: 

1. The stars span a wide range of metallicities between 
[Fe/H] = -0.72 up to super-solar metallicity of [Fe/H] « 
+0.54. 

2. The mean metallicity of the 14 microlensed dwarf and sub- 
giant stars is ([Fe/H]) = -0.08 + 0.47 in good agree- 
ment with the 204 giant stars in Baade's window from 
IZoccali et alj |2008) that have an average metallicity of 
-0.04 dex. However, a two-sided KS-test gives only a low 
30 % probability that microlensed dwarf stars and giant stars 
in Baades window have the same MDFs. The low probabil- 
ity is due to the skewed and uneven metallicity distribution 
of the dwarf stars, with excesses at both low and high metal- 
licities. More observations of microlensed dwarf stars will 
certainly refine the comparison. It is clear thou gh that the 
extreme ly metal-rich MDF for the Bulge that ICohen et alj 
(2009b) propose is not borne out by the larger sample pre- 
sented here. 

3. The abundance trends that the microlensed dwarf stars 
show are surprisingly well-defined. At sub-solar [Fe/H] they 
are more or less coincident with the abundance trends of 
the Galactic thick disc as trac ed by nearby dwarf stars 
dBensbv et alJl2003Ll2005Ll2007h . At super-solar [Fe/H] they 
follow the trends we see for nearby thin disc dwarf stars. 
However, due to the high ages that some of the Bulge stars 
possess at super-solar [Fe/H], and due to the lack of Bulge 
stars at sub-solar [Fe/H] with thin disc abundance ratios we 
see no obvious connection between the Bulge and the thin 
disc. 

4. All stars with sub-solar [Fe/H] are old (around lOGyr) and 
have high [ff/Fe] ratios, consistent with fast enrichment by 
core-collapse supernovae during the early stages of the for- 
mation of the Galaxy. At super-solar [Fe/H] we have a few 
old stars but also three stars with ages lower than 5 Gyr. This 
is inconsistent with, e.g., recent CMDs of field stars in the 
Bulge based on deep imaging with HST/ACS, that show no 
evidence for a young stellar component in the Bulge. The 
average age for our sample of microlensed dwarf stars is 
8.4 + 3.3 Gyr. 

5. Additionally, our results indicate that the red clump stars in 
the Bulge have (V - I) = 1.08. 



Based on these results and conclusions we speculate on the 
origin of the Bulge and we must conclude that it is still poorly 
constrained. 
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Table 5. Elemental abundance ratios, errors in the abundance ratios, and number of lines used, for 13 of the 15 microlensed dwarf 
stars ' . 





[Fe/H] 


[O/Fe]* 


[Na/Fe] 


[Mg/Fe] 


[Al/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Cr/Fe] 


[Ni/Fe] 


[Zn/Fe] 


[Y/Fe] 


[Ba/Fe] 


OGLE-2009-BLG-076S 


-0.72 


0.50 


0.11 


0.36 


0.30 


0.24 


0.34 


0.30 


- 


0.07 


- 


- 


-0.07 


MOA-2009-BLG-493S 


-0.70 


0.43 


0.13 


0.37 


0.28 


0.18 


0.17 


0.34 


0.01 


0.04 


0.25 


- 


-0.07 


MOA-2009-BLG-133S 


-0.64 


0.47 


0.13 


0.39 


0.28 


0.19 


0.21 


0.26 


-0.01 


0.06 


— 


— 


-0.19 


MOA-2009-BLG-475S 


-0.52 


0.33 




0.25 


0.23 


0.18 


0.14 


0.33 




0.07 






-0.18 


MACHO- 1 999-BLG-022S 


-0.49 




0.01 


0.35 


0.31 


0.25 


0.18 


0.25 


-0.03 


0.00 


- 


0.01 


0.15 


OGLE-2008-BLG-209S 


-0.30 


0.34 


0.08 


0.34 


0.24 


0.20 


0.18 


0.29 


0.05 


0.05 


0.16 


-0.08 


0.05 


MOA-2009-BLG-489S 


-0.18 


0.18 


-0.01 


0.24 


0.12 


0.07 


0.09 


0.16 


0.04 


0.01 


0.17 


-0.05 


-0.07 


MOA-2009-BLG-456S 


0.12 


-0.12 


-0.03 


0.09 


0.08 


-0.01 


0.00 


0.08 


-0.06 


-0.03 


-0.03 


_ 


-0.08 


MOA-2008-BLG-311S 


0.36 


-0.27 


0.12 


0.03 


0.08 


0.05 


0.00 


0.04 


-0.01 


0.12 


0.06 


0.11 


-0.15 


MOA-2008-BLG-310S 


0.42 


-0.27 


0.08 


0.08 


0.10 


0.02 


-0.03 


0.09 


0.05 


0.06 


-0.03 


0.06 


0.03 


OGLE-2007-BLG-349S 


0.42 


-0.14 


0.15 


0.05 


0.09 


0.12 


-0.06 


0.02 


-0.03 


0.11 


0.11 


-0.12 


-0.13 


MOA-2008-BLG-099S 


0.44 


-0.10 


0.05 


0.04 


-0.02 


0.00 


-0.05 


-0.01 


-0.16 


-0.02 


0.05 


-0.04 


0.00 


OGLE-2006-BLG-265S 


0.47 


-0.25 


0.14 


0.13 


0.08 


0.06 


-0.05 


-0.06 


-0.07 


0.07 


-0.05 


0.02 


-0.08 




C"[Fe/H] 


0"[0/Fe] 


0"[Na/Fe] 


°"[Mg/Fe] 


°" [Al/Fe] 


C [Si/Fe] 


°"[Ca/Fe] 


°" [Ti/Fe] 


C"[Cr/Fe] 


V [Ni/Fe] 


0"[Zn/Fe] 


C"[Y/Fe] 


C"[Ba/Fe] 


OGLE-2009-BLG-076S 


0.07 


0.17 


0.04 


0.06 


0.05 


0.06 


0.05 


0.03 


- 


0.09 


- 


- 


0.10 


MOA-2009-BLG-493S 


0.14 


0.51 


0.07 


0.10 


0.08 


0.19 


0.30 


0.08 


0.10 


0.17 


0.38 


- 


0.18 


MOA-2009-BLG-133S 


0.17 


0.62 


0.13 


0.21 


0.07 


0.12 


0.36 


0.06 


0.16 


0.15 


— 


— 


0.51 


MOA-2009-BLG-475S 


0.17 


0.39 




0.09 


0.09 


0.10 


0.18 


0.08 




0.15 






0.31 


MACHO- 1 999-BLG-022S 


0.14 




0.07 


0.13 


0.10 


0.13 


0.13 


0.24 


0.14 


0.10 


0.19 


0.31 


0.12 


OGLE-2008-BLG-209S 


0.06 


0.22 


0.07 


0.10 


0.04 


0.08 


0.12 


0.08 


0.09 


0.07 


0.18 


0.20 


0.08 


MOA-2009-BLG-489S 


0.11 


0.25 


0.12 


0.07 


0.05 


0.09 


0.12 


0.08 


0.06 


0.08 


0.09 


0.34 


0.10 


MOA-2009-BLG-456S 


0.09 


0.21 


0.05 


0.08 


0.05 


0.07 


0.10 


0.05 


0.09 


0.08 


0.17 


- 


0.09 


MOA-2008-BLG-311S 


0.07 


0.18 


0.12 


0.09 


0.05 


0.04 


0.11 


0.10 


0.07 


0.05 


0.13 


0.19 


0.11 


MOA-2008-BLG-310S 


0.08 


0.16 


0.10 


0.09 


0.07 


0.05 


0.10 


0.08 


0.06 


0.05 


0.13 


0.23 


0.09 


OGLE-2007-BLG-349S 


0.08 


0.19 


0.15 


0.14 


0.08 


0.07 


0.12 


0.07 


0.05 


0.04 


0.10 


0.17 


0.08 


MOA-2006-BLG-099S 


0.10 


0.20 


0.11 


0.10 


0.04 


0.07 


0.11 


0.11 


0.08 


0.06 


0.10 


0.30 


0.10 


OGLE-2006-BLG-265S 


0.06 


0.17 


0.11 


0.10 


0.05 


0.06 


0.11 


0.07 


0.06 


0.06 


0.43 


0.23 


0.08 




0"[Fe/H] 


°"[0/HJ 


C"[Na/H] 


°"[Mg/H] 


0"[A1/H] 


0"[Si/H] 


C"[Ca/H] 


C"[Ti/H] 


°"[Cr/H] 


f[Ni/H] 


0"[Zn/H] 


0"[Y/H] 


0"[Ba/H[ 


OGLE-2009-BLG-076S 


0.07 


0.10 


0.05 


0.06 


0.05 


0.05 


0.10 


0.08 


- 


0.06 


- 


- 


0.11 


MOA-2009-BLG-493S 


0.14 


0.38 


0.09 


0.16 


0.07 


0.06 


0.43 


0.13 


0.07 


0.05 


0.26 


- 


0.18 


MOA-2009-BLG-133S 


0.17 


0.46 


0.05 


0.36 


0.18 


0.06 


0.52 


0.12 


0.07 


0.04 


— 


— 


0.36 


MOA-2009-BLG-475S 


0.17 


0.23 




0.17 


0.13 


0.08 


0.31 


0.14 




0.07 






0.21 


MACHO-1999-BLG-022S 


0.14 




0.14 


0.23 


0.06 


0.06 


0.25 


0.16 


0.11 


0.08 


0.10 


0.22 


0.15 


OGLE-2008-BLG-209S 


0.06 


0.16 


0.11 


0.15 


0.07 


0.03 


0.17 


0.06 


0.06 


0.03 


0.16 


0.16 


0.08 


MOA-2009-BLG-489S 


0.11 


0.15 


0.13 


0.14 


0.07 


0.03 


0.21 


0.09 


0.07 


0.04 


0.06 


0.29 


0.09 


MOA-2009-BLG-456S 


0.09 


0.14 


0.07 


0.12 


0.09 


0.04 


0.15 


0.10 


0.14 


0.04 


0.15 


_ 


0.10 


MOA-2008-BLG-311S 


0.07 


0.12 


0.17 


0.14 


0.10 


0.06 


0.17 


0.07 


0.06 


0.04 


0.07 


0.13 


0.07 


MOA-2008-BLG-310S 


0.08 


0.10 


0.14 


0.13 


0.09 


0.04 


0.15 


0.08 


0.05 


0.05 


0.08 


0.21 


0.12 


OGLE-2007-BLG-349S 


0.08 


0.13 


0.19 


0.18 


0.08 


0.04 


0.17 


0.09 


0.07 


0.06 


0.10 


0.17 


0.12 


MOA-2006-BLG-099S 


0.10 


0.11 


0.18 


0.18 


0.08 


0.05 


0.20 


0.07 


0.06 


0.07 


0.09 


0.25 


0.09 


OGLE-2006-BLG-265S 


0.06 


0.12 


0.14 


0.14 


0.07 


0.04 


0.14 


0.08 


0.07 


0.04 


0.43 


0.21 


0.08 






N 


^Na 


NfAg 


N A \ 


Wsi 


Afca 


Af Ti 


Na 


^Ni 


N Zn 


N Y 


M5a 


OGLE-2009-BLG-076S 


57 


3 


1 


5 


4 


15 


10 


2 


- 


11 


- 


- 


2 


MOA-2009-BLG-493S 


80 


3 


1 


6 


6 


19 


10 


5 


2 


21 


1 


- 


4 


MOA-2009-BLG-133S 


68 


2 


1 


4 


6 


17 


10 


4 


1 


15 


— 


— 


3 


MOA-2009-BLG-475S 


53 


1 




5 


3 


12 


9 


2 




14 






3 


MACHO-1999-BLG-022S 


97 




3 


1 


2 


7 


16 


14 


2 


21 




2 


3 


OGLE-2008-BLG-209S 


146 


3 


4 


5 


7 


27 


19 


36 


10 


40 


3 


4 


4 


MOA-2009-BLG-489S 


114 


3 


2 


6 


6 


27 


13 


16 


6 


37 


3 


2 


4 


MOA-2009-BLG-456S 


91 


3 


2 


5 


6 


26 


11 


7 


3 


33 


3 




2 


MOA-2008-BLG-311S 


118 


3 


4 


5 


5 


26 


17 


14 


9 


36 


1 


2 


3 


MOA-2008-BLG-310S 


122 


3 


4 


4 


7 


27 


17 


22 


13 


42 


1 


4 


3 


OGLE-2007-BLG-349S 


103 


3 


4 


3 


3 


23 


18 


24 


9 


39 


3 


5 


4 


MOA-2006-BLG-099S 


119 


3 


4 


4 


4 


25 


18 


30 


12 


37 


3 


5 


4 


OGLE-2006-BLG-265S 


92 


3 


4 


2 


3 


23 


16 


13 


8 


30 


2 


2 


3 



' No abundances are given for MOA-2009-BLG-259S as the errors in the stellar parameters were too large , and for OGLE-2007-BLG-514S we only redeteremined 

stellar parameters and [Fe/H]. Abundance ratios for OGLE-2007-BLG-514S can be found in Epstein et al. (2009). 

* Note that the abundance ratios for oxygen have been corrected for NLTE effects according to the empirical formula given in iBensbv et ail 120041) . 



